Description beyond the mean field approximation of an electrolyte confined between 

two planar metallic electrodes 
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We study an electrolyte confined in a slab of width W composed of two grounded metallic parallel 
electrodes. We develop a description of this system in a low coupling regime beyond the mean field 
(Poisson-Boltzmann) approximation. There are two ways to model the metallic boundaries: as ideal 
conductors in which the electric potential is zero and it does not fluctuate, or as good conductors in 
which the average electric potential is zero but the thermal fluctuations of the potential are not zero. 
This latter model is more realistic. For the ideal conductor model we find that the disjoining pressure 
is positive behaves as 1/W 3 for large separations with a prefactor that is universal, i.e. independent 
of the microscopic constitution of the system. For the good conductor boundaries the disjoining 
pressure is negative and it has an exponential decay for large W . We also compute the density 
and electric potential profiles inside the electrolyte. These are the same in both models. If the 
electrolyte is charge asymmetric we find that the system is not locally neutral and that a non-zero 
potential difference builds up between any electrode and the interior of the system although both 
electrodes are grounded. 

PACS numbers: 61.20.Qg, 82.45.Gj, 82.45.Fk 
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I. INTRODUCTION 

In this paper we study an electrolyte solution confined 
between two parallel planar metallic electrodes. The 
study of the electrical double layer near an electrode and 
more generally near any object submerged in an elec- 
trolyte is of crucial importance in chemical physics and 
in colloidal science. This problem was first considered 
by Gouy [l| and independently by Chapmann 2] almost 
a century ago. Their work is part of the foundations of 
colloidal science Q and the physics of electrolytes 0. 

However their work and its developments are based on 
a mean field description: the Poisson-Boltzmann equa- 
tion. Although this mean field approach describes accu- 
rately several properties of the systems, in some situa- 
tions it misses some subtle effects due to correlations. As 
an example we can mention the old controversy about 
the possibility of attraction between charged-like col- 
loids y, la] recently renewed by some experimental re- 
sults ja, H It has been shown 0, 0, El that 
the mean field approach (actually any local density ap- 
proximation) cannot predict any attractive effective in- 
teraction. Therefore the study of electrolyte suspensions 
beyond the mean field approximation is important. 

This paper is oriented in that sense, although we will 
not consider here the problem of charge-like attraction 
between colloids, but the study of an electrolyte solution 
confined between two parallel metallic planar electrodes 
beyond the mean field approximation. We will be inter- 
ested in questions like what is the force exerted on the 
planar electrodes, it is attractive or repulsive, etc. . . ? To 
have a clear picture of the role of the correlations in this 
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problem we will consider the case when the two electrodes 
are grounded. The mean field picture in this case is very 
simple: the mean field potential in the electrolyte is zero 
everywhere and the fluid is uniform and locally neutral. 
We will describe the first fluctuations around this mean 
field picture in a low coupling regime where the average 
coulombic energy of the microions of the solution is much 
smaller than their thermal energy. 

We should mention that this same problem was re- 
centlyconsidered in Ref. However the authors of 

Ref. [13j made a mistake that has lead them to the wrong 
conclusions, as explained below. The electrolyte is con- 
fined between two conductor parallel planes. Each par- 
ticle polarizes the planes. There is an interaction energy 
between each particle and the polarization charge that it 
induces in the electrodes. In Ref. the authors forgot 
to include this energy in the hamiltonian and this error 
make most of their conclusions incorrect. 

The outline of this paper and our main results can be 
summarized as follows. In Sec. [H] we present the models 
under consideration and explain the method 0] used to 
find the thermodynamic properties of the system. Actu- 
ally there are two ways to model the metallic electrodes. 
In a first model the boundaries are supposed to be made 
of an ideal conductor material. This model is very simple 
but it has the defect of neglecting the fluctuations of the 
electric potential in the electrodes. In the second model, 
which is more realistic, the electrodes are supposed to 
be genuine Coulomb systems with very good conducting 
properties. We will call this model the good conductor 
model. In this model the screening length vanishes inside 
the conducting electrodes, the average electric potential 
is zero in the electrodes, but the electric potential can 
fluctuate. We will obtain results for both models and 
compare them using results from Ref. |15|. 

In Sec. IIIII we compute the grand potential of the sys- 
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tern and the pressure. For the ideal conductor model we 
find that the disjoining pressure is positive. For large sep- 
arations W of the slab, the disjoining pressure behaves 
as 1/W 3 . On the other hand, for the good conductor 
model, the disjoining pressure turns out to be negative 
and it decays as e~ 2KW for large separations, with k the 
inverse Debye length. 

Finally in Sec. II VI we find the microion density profiles 
and the electric potential inside the electrolyte. These 
quantities are the same regardless of the model used to 
describe the electrodes (ideal conductor or good con- 
ductor). In that section we retrieve the important re- 
sult ee mm 

that for charge asymmetric electrolytes, 
a non-zero potential difference builds up between each 
electrode and the middle of the electrolyte solution and 
the system is not locally neutral although both confining 
plates are grounded. 

II. MODEL 

As explained in the Introduction, the system un- 
der consideration is an electrolyte confined between two 
grounded conductor planar electrodes separated by a dis- 
tance W. Let us choose the x-axis in the direction per- 
pendicular to the electrodes, the origin is in the mid- 
dle of the electrodes and the electrodes are located at 
x = ±W/2. We will eventually also consider the limiting 
case when W — + oo. In this case we shall use the coor- 
dinate X = x + W/2 which measures the distance from 
one electrode. 

We will model the electrolyte as been composed of sev- 
eral species of point-like microions with charges q a la- 
beled by a Greek index. The solvent will be modeled as 
a continuum medium with dielectric constant e. As it is 
well-known [ltA |20| a system of charged point particles 
described by classical statistical mechanics is not stable. 
In principle we should introduce some short-distance cut- 
off, for instance the radius of the ions, to avoid the col- 
lapse of particles of different sign. However in the low 
coupling regime considered here most physical quantities 
have a well defined value when this short-distance cutoff 
vanishes. Therefore our model will give valuable quanti- 



tative information on the system provided the ion radius 
is much smaller that the others lengths involved in the 
problem. 

The position of the z-th particle of the species a will 
be labeled as r Q j. We shall work in the grand-canonical 
ensemble at a reduced inverse temperature (3=1/ (fc^T), 
with ks the Boltzmann constant and T the absolute tem- 
perature. The average number of particles (N a ) of the 
species a is controlled by the chemical potential fi a . We 
shall use the fugacity £ Q = e l3Ua /A^ where A a is the 
thermal de Broglie wavelength of the particles which ap- 
pears as usual in classical (i.e. non-quantum) statistical 
mechanics after the trivial Gaussian integration over the 
kinetic part of the hamiltonian. We shall impose the 
pseudo-neutrality condition 

]>>«C« = 0. (2.1) 

a 

In the appendix B of Ref. E3 it is explained that this 
choice is equivalent to suppose that there is no electric 
potential difference between the plates and the interior 
of the system in the mean field approximation. The 
pseudo-neutrality condition ensures the neutrality of the 
reservoir, however as we will see below it does not en- 
sure the global neutrality of the system confined with 
Dirichlct boundary conditions. As we will show in sec- 
tion ]W\ there is a non- vanishing charge density near the 
electrodes. However, this charge density induces in the 
electrodes a polarization charge of opposite sign and the 
total charge (system plus electrodes) is zero. 

There are two models that can be used to describe the 
electrodes |2l|. The simpler one is the ideal con- 
ductor model. In this model, the interaction potential 
between two unit charges located at r — (x, y, z) and 
r' = (x' , y', z') is the solution of Poisson equation 

Av(r,r') = -^-5(r-r') (2.2) 

satisfying the Dirichlet boundary conditions t>(r, r') = 
if x' = ±W/2. It can be computed using, for example, 
the method of images, 



v(r, r') 



[(x -x' + 2nW) 2 + (r ± - r' ± ) 2 ] 1/2 [(x + x' + (2n + l)W) 2 + (r± - r' ± ) 2 ] 1/2 



(2.3) 



with rj_ = (y, z) the transversal part of the position vec- unconfined system 
tor r and e is the dielectric constant of the solvent. For 

future reference we define the Coulomb potential for an n . 1 1 . . 

u u (r,r') = — j (2.4) 



e r — : 



which will be needed in the following. 
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Notice that with this ideal conductor model, the mi- 
croscopic electric potential deduced from (|2.3|) vanishes 
on the electrodes, for any configuration of the system. 
Then, not only the average electric potential will vanish 
on the electrodes, but also its fluctuations. 

In a real situation, even if the electrodes have very 
good conducting properties, the electric potential inside 
the conductor will be subjected to thermal fluctuations. 
This leads us to the other model that can be used to 
describe the electrodes, the good conductor model. In 
this model the electrodes are genuine Coulomb systems 
with good conducting properties. We will restrict our- 
selves to a classical (non-quantum) description of these 
conductors. We will consider that the screening length in 
the conductors vanishes. This ensures that the electrodes 
have good conducting properties. It is shown in |l5j that 
in this limit the average electric potential vanishes on the 
electrodes. However the fluctuations of the potential do 
not vanish in this limit. 

Then a natural question arises: how are related the 
ideal conductor model and the good conductor one? This 
question has already been addressed in Ref. |l5j and we 
will make extensive use of the results presented there. It 
is shown in Ref. that the densities and electric po- 
tential profiles inside the electrolyte are the same in both 
models. Thus we will work with the (more simpler) ideal 
conductor model for the determination of these profiles 
in Sec. |TV| 

On the other hand, the results for the grand poten- 
tial and for the pressure are different in both models, 
but there are simple relations between them [l5|. We 
will use the ideal conductor model for the calculation 
of these quantities in Sec. IHII and using the relations 
found in [15j we will obtain the results for good conduc- 
tor model. From now on we consider the ideal conductor 
model unless stated otherwise. 

Although to write down the hamiltonian of the system 
is a trivial exercise in electrostatics, to clearly show what 
the problem is with the previous study |13| of this system, 
we will detail a few (well-known) points before proceed- 
ing. First, consider the case when only a planar electrode 
is located at X — 0. Bringing first a unit charge from 
infinity to a position r = (X, y, z) at a distance X from 
the plane cost a non-zero energy, contrary to the case 
of an unconfined system. This is because of the interac- 
tion between the particle and the polarization charge it 
induces in the plane. In this very simple geometry this 
interaction can also be understood as the potential en- 
ergy between the particle and an image charge located at 
r* = (— X,y, z). This energy is — l/{AeX) which can be 
formally written as (l/2)[w(r, r) — w°(r,r)] (in this case 
v(r, r') is the potential (|r — r'| _1 — |r — r* \~ l )/e when 
only one electrode is present). This interaction energy 
should be included in the hamiltonian. 

Following the same lines, in the general case of two 



metallic planes the potential energy of the system reads 

a, 7 

+ [v(r a ,i,r a ,i) -v°(r aji ,r aii )] . 

a i 

In the first sum the prime means that the case a = 7 
and i = j should be omitted. The second sum is the 
energy between each particle and the polarization charge 
it has induced in the electrodes as discussed previously. 
Introducing the microscopic charge density defined as 

p(r)=^^ gtl 5(r-r a , I ) (2.6) 

a i 

we can formally write the potential part of the Hamilto- 
nian of the system as 

H = i J dr J dr' p(t)v(t,t')p(t') (2.7) 

a i—1 

The domain of integration in the first term is the space 
between the two parallel electrodes (-W/2 < x < W/2). 
Notice that from the first term written in terms of "con- 
tinuous" fields we subtract the infinite "self-energy" of a 
particle i>°(r, r) but with the potential energy v° corre- 
sponding to an unconfined system. 

With the Coulomb potential v given by Eq. (|2.4H this 
self-energy is infinite. The appearance of infinite quanti- 
ties can be avoided and the subsequent analysis can be 
done more rigorously by modifying the Coulomb poten- 
tial v°(r, r') by introducing a short-distance cutoff which 
smears out the singularity at r = r'. Then, at the end of 
the calculations one can take the limit of a vanishing cut- 
off. This has been done rigorously in Ref. [2^J . It turns 
out that the final results are the same as if one takes the 
short-distance cutoff equal to zero from the start. There- 
fore to simplify the algebra we will take from the start a 
system of point particles without short-distance cutoff. 

Now we follow the method proposed recently by the 
author and collaborators in Refs. [Fa . 1231 ] to study in gen- 
eral confined Coulomb systems in a low coupling regime. 
Let us define the coulombic couplings r a = /3q%,Ca / 'e- 
The method proposed in Ref. [Tj| is valid for r a C 1. 

In the method exposed in Ref. 0] the sine-Gordon 
transformation |24| is performed in the grand-canonical 
partition function, then the action of the corresponding 
field theory is expanded to the quadratic order (valid in 
the low coulombic coupling regime) around the stationary 
(mean field) solution (here = 0). For details the reader 
is referred to Ref. 01 an d to the appendix Then the 
grand partition function can be written as 




V<t> exp[-5(0)] (2.8) 
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with 



T>4> exp 



(r) 



(3eA 



4tt 



(r) dr 



(2.9) 



and S(4>) is an action, quadratic in <j>, given by 



(2.10) 

where :: 4>(y) 2 :: is a pseudo-normal ordered product, see 
the appendix 1X1 for details. 

The field <f>{r) is a mathematical intermediary. At the 
mean field level, the stationary equation for the action 
(before it is expanded to the quadratic order) is Poisson- 
Boltzmann equation, and i(f>(r) can be interpreted as 
the electric potential, however this relation breaks down 
when we consider the fluctuations as in the present case, 
for instance the correlations of <p(r) are short-ranged [2f| 
whereas the correlations of the electric potential are 
known to be long-ranged j2|| [2?]]. The Gaussian func- 
tional integration in Eq. (|2.8fl can be performed to 
obtain 



n -f n 



-1/2 



(2.11) 



where A„ are the eigenvalues of the Laplacian operator 
satisfying the Dirichlet boundary conditions and A^ are 
the eigenvalues of the Laplacian operator defined in the 



whole space R 3 without boundaries. We will call this case 
in the following the free boundary conditions case. The 
volume of the system is V and K = \/J2a ^ n Ca(3q 2 /e 
is the inverse Debye length. The second product in 
Eq. p.llfl involving comes from the subtraction of 
the self-energy term v°(r, r). 



III. GRAND POTENTIAL AND PRESSURE 

A. Ideal conductor model 

In this part we use the ideal conductor model to de- 
scribe the electrode. As a remainder of this, we will use 
a superscript "id" in the thermodynamic quantities. 



1. Grand potential 

For the present geometry the eigenvalues of the Lapla- 
cian for Dirichlet boundary conditions and free boundary 
conditions respectively are A = — k 2 — (mr) 2 /W 2 with 
neN* and k e R 2 and A°„ = -K 2 with K € R 3 . 
We find that the grand potential f2 ld takes the form 
Q id = Q ig + nl d xc with n ig = k B TVJ2 a C« the ideal gas 
contribution and the excess grand potential. From 
Eq. H2.11fl we find the excess grand potential per 
unit area of a plate 



J 



0w, 



id 



2(2tt) 



mn 



d 2 k- 



The product under the logarithm can be performed exactly |23| to obtain 



In 



k sinh(WVn 2 + k 2 ) 



V* 2 



sinh(fcVF) 
I 



W K 

2(2tt) 



k dk 



Wk 2 



d 3 K 
K 2 



dK . 



(3.1) 



(3.2) 



Notice that we introduced two ultraviolet cutoffs /c max 
and K m &x for both integrals since each integral, taken 
separately, is ultraviolet divergent. However together 
they should give a finite result when fc max — > oo and 
-Kmax - * oo as far as the bulk properties are concerned. 
Indeed, in the limit W — *■ oo we should recover the well- 



known bulk result Hi H2 f3uj b = -k 3 VT/(12tt). This 
requirement imposes that the cutoffs should be related 
by i^max — Then doing the change of vari- 

able K = irk/2 in the second integral the excess grand 
potential per unit area can finally be written as 



J 



id 

exc 



1 

47T 



k In 



i1i(WVk 2 + k 2 ) 



V^ 2 



sinh(fcVF) 



n 2 W 



dk. 



(3.3) 



In principle we should take the limit fc max — + oo, how- 
ever it should be noted that the above expression has a 



logarithmic divergence when fc n 



oo which manifests 
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itself in the surface tension. This can be seen clearly if 
we expand w exc for kW 3> 1, 



#4 



d 

exc 



K 3 W 



2/? 7 - 



C(3) 



12tt ' ' ' 16nW 2 
with the surface tension 7 given by 

,2 / . j 



0(e 



-2kW\ 



Pi 



167T V U 



(3.4) 



(3.5) 



and C(3) is the Riemann zeta function evaluated at 3 (not 
to be confused with the fugacities). In Eq. i|3.4|) all terms 
that vanish when fc max — *■ 00 have been omitted. 

A few comments are in order. Concerning the surface 
tension 7 it is divergent when the cutoff k max — » 00. This 
is normal: it is due to the strong attraction that each 
particle and its images of opposite charge in the elec- 
trodes feel. The small coupling regime of an electrolyte 
near a plane metallic wall can also by studied from a dia- 
grammatic Mayer expansion. This is done in section 5 of 
Ref. [2£j for a two-dimensional Coulomb system. These 
calculations can easily be adapted to a three dimensional 
system to show that the surface tension 7 is related to 
the integral of the screened interaction energy between a 
particle and its image: — exp(— 2kX)/(4X). This energy 
is not integrable at short distances and its integral has a 
logarithmic divergence at X = 0. In this picture one can 
impose a short-distance cutoff D: the particles cannot 
approach below this distance to the electrode, then the 
surface tension is proportional to hinD. Actually our 
ultraviolet cutoff k max ocl/D. 

The second comment concerns the algebraic finite- 
size correction fcsT£(3) / (WnW 2 ) to the grand potential. 
This finite-size correction is universal, it does not depend 
on the details of the microscopic constitution of the sys- 
tem, and it has been proved to exist even beyond the low 
coupling regime considered here provided that the elec- 
trolyte is in a conducting phase and it has good screening 
properties, in particular if it can screen an external in- 
finitesimal dipole [3(]] ■ This correction is a consequence 
only of the screening properties of the system, that ex- 
plains its universality. However as we will see later, this 



term is only present for the ideal conductor model con- 
sidered in this section. For the good conductor model 
this algebraic correction in the grand potential is absent. 

We should also mention that this algebraic finite- 
size correction is not present in the case of insulating 
plates j33, El E3 • Indeed if the boundary is made of a 
dielectric material there is a subtle cancellation between 
the term found here and the one from the Lifshitz inter- 
action (the Casimir effect). For dielectric boundary con- 
ditions the disjoining pressure has an exponential de cay 
6 -2kW a |- j ar g e separations W and it is attractive [3lll32l| 
as opposed to the present case of ideal conductor bound- 
ary conditions. Further details on this interesting differ- 
ence are discussed in the appendix FBI 
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FIG. 1: The disjoining pressure of the electrolyte confined 
by ideal conductor electrodes. It is positive and always de- 
creasing with increasing W indicating that there is a repulsive 
force between the two ideal conductor parallel plates. 



2. Pressure 



The pressure is obtained from the usual relation p = 
duj/dW. From Eq. (|3.3|) we find that the excess pres- 
is given by 



sure p^c 



J 



Pp, 



id 

exc 



4?r ,,o 



r« 2 



k 2 coth(fcVF) - k\Jk 2 + k 2 coth(wVfc 2 + k 2 ) 



r 



dk. 



(3.6) 



Although the grand potential has an ultraviolet diver- 
gence and should be regularized as explained earlier, the 
pressure proves to be well defined for /c max — > 00 (and 
W 7^ 0). This is expected since from the large- W ex- 
pansion (|3.4|) of the grand potential we can see that the 
ultraviolet divergent part (the surface tension contribu- 
tion) docs not depend on W . Notice however that for 



W — ► the pressure is divergent. Let us mention that the 
non-divergence of the pressure with the cutoff and more 
precisely the fact that it is independent of the surface 
tension 7 is special to this planar geometry. If we were 
to consider a confining geometry with curved boundaries 
(for example an electrolyte confined in a spherical do- 
main) the surface tension would be a dominant term in 
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the pressure: due to the curvature R the disjoining pres- 
sure for large systems would be pd oc —7 / R, see Ref. 
for an example of this effect. 

Doing a few manipulations to Eq. I|3.6[) we can cast the 
pressure in a form more adequate to study the disjoining 
pressure p^ d , difference between the pressure p ld and the 
bulk pressure p b , and its large- W behavior. The bulk 
pressure, expressed in terms of the fugacities, is obtained 
from the limit W — > 00 of Eq. I|3.4[) . and it is given by 



The well-known expression of the bulk pressure in terms 
of the densities will be recovered in the next section, 
Eq. 1)4. f 7|l . when we obtain the expression of the bulk 
densities in terms of the fugacities. 



0p b = ^2ta 



K 

12tt 



(3.7) 



Then we find the disjoining pressure 



id 



C(3) 
8irW 3 

C(3) 
8ttW 3 



K 

47T 



i\Ju 2 + 1 



1 — coth 



(kW V~i 



du 



+ 0(e- 2KW ) 



(3.8) 
(3.9) 



For large separations W of the electrode plates the dis- 
joining pressure is positive and decays as l/W 3 . The 
force exerted on the plates is therefore repulsive and 
the system is stable: if it is compressed the pressure in- 
creases. This actually holds for any separation as it can 
be seen in Fig. ^ 



B. Good conductor model 

We now study the electrolyte confined by good con- 
ductor electrodes. We will not use any superscript in the 
thermodynamic quantities in order to differentiate them 
from the ones computed for the ideal conductor model. 



1. Free energy and grand potential 

Let F be the free energy of the system composed by 
the electrolyte and the electrodes. These later ones are 
considered as Coulomb systems with a vanishing screen- 
ing length (high density of charged particles) in the region 
x < —W/2 and x > W/2. The electrolyte is in the region 
-W/2 <x< W/2. 

For the following discussion it is useful to introduce 
Fq the free energy of the electrodes in absence of the 
electrolyte. It is shown in Ref. [H| that F = F id + F 0> 
where F ld is the free energy of the electrolyte in the ideal 
conductor model, which can be obtained from the results 
of the previous section. The argument of Ref. can 
easily be applied to the grand potential, thus 



LO 



(3.10) 



with uj the grand potential (per unit area of the plates) 
of the system when the electrodes are described as good 



conductors, a-> ld the grand potential (per unit area) for 
the ideal conductor model, which can be obtained from 
Eq. (|3.3I) . and c^o the grand potential of the good con- 
ducting electrodes alone in the space in the absence of 
the electrolyte. 

The large- W expansion of this latter term loq is given 
by Lifshitz theory |34[ in the classical regime, it is 
— ksT^S) / (16ttW 2 ) . Thus it cancels the similar contri- 
bution that we found in uj ld in the last section, Eq. (|3.4|l . 
When the electrodes are described as good conductors 
(which is a more realistic model) there is no long-range 
contribution in 1/W 2 to the grand potential of the sys- 
tem as opposed to the ideal conductor model. 



2. Pressure 

The pressure is the force per unit area that the elec- 
trolyte exerts on one plate, say the one at x = W/2. It 
can be computed by means of the Maxwell stress tensor 
Tpv. It is p = —T xx evaluated at x = W/2. In Ref. [Ht 
it is shown that the stress tensor in the ideal conductor 
model Tji and the one in the good conductor model T^ v 
are related by 



-in 



dud' 



(3.11) 

with G*(r,r') = v°(r, r') — v(r,r') is (minus) the "im- 
ages" contribution to the Coulomb potential v in the ideal 
conductor model. 

The r.h.s. of Eq. (|3.11J) can be computed explicitly giv- 
ing 



J-xx ~ 1 xx 



C(3) 
8irW 3 



(3.12) 
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FIG. 2: The disjoining pressure of the electrolyte confined by 
good conductor electrodes. It is negative and always increas- 
ing with increasing W indicating that there is an attractive 
force between the two conductor parallel plates. 



Thus the pressure for the good conductor model is 

C(3) 



P = P 



8nW 3 ' 



(3.13) 



Using Eq. I|3.8|l we finally obtain the disjoining pressure 
when the electrodes are modeled as good conductors 



1 



du . 



K 3 i-co 

(3pd = — / u\J u 2 + 1 1 - coth ( kW \/i 
4tt J I V 

(3-14) 

Notice that since the function coth in the integrand is 
greater than 1, the disjoining pressure is always negative. 



Figure[5]shown a plot of the disjoining pressure as a func- 
tion of the width W. We notice that the pressure is now 
an increasing function of W. This is just the opposite 
behavior that the one obtained with the ideal conductor 
model. Also for large-W the disjoining pressure decays 
exponentially as e" 



-2kW 



IV. DENSITY AND ELECTRIC POTENTIAL 
PROFILES 



The difference in the results for the pressure using the 
ideal conductor model and the good conductor one are 
drastic. However the results for the density and electric 
potential profiles inside the electrolyte are the same in 
both models as it was shown in Ref. |ll|. Therefore we 
will concentrate in this section on the ideal conductor 
model which is more tractable. 



A. Density 

The density n a (r) can be obtained from the usual func- 
tional derivative 



n a {r) = ( a (r) 



ShxE 



In the appendix 1X1 it is shown that 



(4.1) 



J 



n a (r) = ( a 



l-^[VDH(r,r) 



(3 2 q a y-, 

7 



v°(r, v))+^Y. Cr«? / Wr', r) [« DH (r', r') - v°(r', r')] dr> 



(4.2) 



r 



where vj^n(r,r') is the Debye-Huckel potential, solution 
of Debye-Hiickel equation 



A 



47T 



WDH(r,r') = <5(r-r') 



(4.3) 



satisfying the Dirichlet boundary conditions vuu{r, r') = 
if x' = ±W/2. Eq. (|4.2|) gives the density up to the 



J 



3/2 

order TJ in the coulombic couplings. For the present 
calculations we found that the most convenient form for 
vdh is as a Fourier transform in the transverse direc- 
tion = (y,z). In Fourier transform, Debye-Hiickel 
equation (|4.3II reduces to an ordinary linear differential 
equation in the x variable, which can be easily solved. 
Then we find 



UDH(r,r') = 



4tt f dk sinh (Vk 2 + k 2 (f - x')) sinh 



>w 

, 2 



ik r I 



(2n) 2 



Vk 2 + k 2 sinh (VFa/P 



(4.4) 
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if x < x' and exchange the roles of x and x' if x' < x. Using this expression into (|4.2(l we find that the density can 
be expressed as 



1 



2s 



-fi(nx) + 



f 2 {nx) 



with 



h(Z) = - 



2k 



sinh ( \Jk 2 + 1 



sinh ( \Jk 2 + 1 



sinh [VW+IW 



dk 



and 



/ 2 (i) 



cosh: 



cosh(W72) Jo 



l-^^^cothf^y^TT 



Ak 2 + 3 



A; cosh (2x\Jk 2 + l) 



o VPTT (4/c 2 + 3) sinh (wVFTI 
/ - k coth f W^Vfc 2 + l) - 1 



dk 



dk 



(4.5) 

(4.6) 

(4.7) 
(4.8) 
(4.9) 



where we have used distances measured in Debye length units x — kx and W = kW. Notice that the factors 

TT^T, 3/2 

multiplying /j and 72 in the density, Eq. (|4.5(l . are of order T Q in the coulombic couplings. Our approach neglects 

3 

corrections of higher order than r a ' 



3/2 



After doing the change of variable u — vP+T in the above integrals some of them can be performed explicitly 
and doing some manipulations we find the following convenient expressions for fi(x) and f2(x) 



(W-2x) ~(W+2x) 

h(x) = 1 + + + 2 

W~2x W + 2x 



00 e -suw cos h(2ui) du 1 



1-e 



-2uW 



W ' 



and 



with 



and 



Hx) = f?\x) + ff\x) - 1 - ± ln(l - e- 2W ) 



f (i) 



cosh: 



cosh(W/2) 



/ 2 (2) (x) 



T-*Ei 



e^+ 5 Ei 



00 „-3uW 



In 3 



'w 



1 (4u 2 -l)(l-e 

- e -(#-*> Ei 
-e-(* +£ ) Ei [- 



( K 



■ G?1i 



(4.10) 
(4.11) 

(4.12) 

(4.13a) 
(4.13b) 
(4.13c) 

li (4u 2 -l)(l-e- 2MW/ ) 

where Ei(z) = — J' x ' z e^ t /tdt is the exponential integral function. The advantage of these latter expressions is that 

one can immediately see that the terms written as integrals are of order (D(e~ 2W ) when W — * 00. Therefore we can 
easily obtain the expression for density in the case of one electrode alone, with X = x + W/2, 



cosh(2ux) du 



n a (X) = ( a 



2e 



,-2kX 



2kX 



In 3 



e KX Ei(-ZkX) - e~ KX ~Ei{-nX) 



(4.14) 
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Far away from the metallic wall, X — > oo, we find the bulk density 



^[ L + 2e e*K 



(4.15) 



Replacing back into Eq. 14.1411 we find an expression for the density profile in terms of the bulk density 



n a {X) 



1 + 



f3qle 



-2k dh X 



AeX 



3„b r 



+ 



ln3\ e KDH * Ei(-3K D H^) - e- KDHX Ei(-K D H^0 



(4.16a) 



(4.16b) 



r 



3/2 

with corrections of higher order than r Q . Here kdh = 

^JAtt[3 n^qZ/e. We recover the expression that Aqua 

and Cornu have previously obtained in their studies of 
the prop erties of a classical Coulomb system near a 
wall llTJ, [18| using diagrammatic methods, u p to a 
small difference: the first term (I4.16a|l appears in 18] as 
an exponential Boltzmann factor of the screened inter- 
action between a particle and its image. Here, this ex- 
ponential Boltzmann factors appears linearized since the 
coulombic coupling is small. Our density profile and the 
one found in |lq will agree, in the low coupling regime 
considered here, for distances not to close to the electrode 
and our results are reliable in this case. On the other 
hand, very close to the electrode, at distances compara- 
ble to the ion radius, our results will differ from those of 
Ref. 0| , this is a defect of the point-like model for the 
microions used here. 

We can use Eq. I|4.15(l which relates the fugacities to 
the bulk densities into the expression l|3.7[) of the bulk 
pressure expressed in terms of the fugacities to recover 
the well-known equation of state of Debye-Hiickel the- 
ory (3 



Returning to the general case, for any arbitrary sepa- 
ration W of the plates it can be noticed that the density 
diverges at x = ±W/2 as l/(x^fW/2). The density does 
not have a finite value at the contact of the electrodes but 
it diverges. This is a expected behavior, since each par- 
ticle is strongly attracted to its images in the electrodes. 
This is related to the divergence of the surface tension 
and the necessity to impose a short-distance minimum 
distance of approach of the particles to the planar elec- 
trodes D cx 1/fcmax as explained in the previous section. 
The logarithmic divergence in In nD of the surface ten- 
sion is closely related to the divergence of the densities 
as l/(x =F W/2) at the contact of each electrode. 

The charge density turns out to be 



p(x) = y^^q a n a (x) 




P(kx) (4.18) 



DH 

24tt 



(4.17) with the reduced charge density 
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p(x) = h(x)+f 2 (x) 

e -(W-2x) e -(W+2x) 



W-2x W + 2x 



i (4u 2 - 1)(1 - e - 2uW ) cosh(W/2) 
cosh(2uai) du 



(4u 2 -l)(l-e- 2 " w ) 



(4.19a) 




(4.19b) 



r 



In the case of a two-component symmetric electrolyte, 
qi = —q2, and we have C-yl^ — 0, therefore the system 
is locally neutral p{x) — 0. For a general asymmetric 
electrolyte Yl a 9oCq and the system is not locally 
neutral. Furthermore the charge density diverges near 
the plates as l/(x=F W/2) which is not integrable. Then 
the total charge induced in the electrodes is infinite if the 
particles are allowed to approach the electrodes as near 
as they can. 

Fig. EH show several charge density profiles for different 
values of W with k fixed. As expected if kW 3> 1 the 
profiles for different values of W are very similar since 
the corrections to the case W — > oo are of order e~ 2nW . 
This can be seen in the plots for kW — 5 and kW = 10 in 
Fig. |3| The differences from the case W — ► oo can be only 
be noticed for small values of kW as in the cases kW = 1 
and kW = 0.16 of Fig- EI However let us remark that for 
any value of W the charge density from an electrode up to 
the middle of the slab is strictly monotonous (increasing 



J 



or decreasing depending on the sign of J2 a (a<la)- 

B. Electric potential 

For the present geometry, the electric potential can be 
computed from the charge density as 

4tt f x 

§(x) - $(0) = — / (x 1 -x)p(x')dx' . (4.20) 
£ Jo 

This gives 



27T/3 

*(a;) - $ = -^f- 



($(kx) - $(0)) (4.21) 



with the reduced electric potential 



- $(0) 



e w72 Ei(-3W72) - er w ' 2 Ei(-iy/2) 
/ W 



-(¥+*) Ei 



e -(¥-x) E i - 




e T-+ £ Ei I -3 ( h / 



W 



In 3 



oo „~3uW 



u 2 e~ 2uW du 



i {4u 2 - 1){1 ~ e~ 2 ^ w ) 



1 — cosh x 
cosh(FF/2) 



(cosh(2ui) — 1) 
i (4u 2 -l)(l-e- 2 ^) 
I 



du . 



(4.22a) 
(4.22b) 

(4.22c) 

(4.22d) 

(4.22e) 



Fig. 0| shows the electric potential profile for different An interesting quantity is the potential difference be- 

values of the width W. tween a plate (for example x — W/2) and the middle of 
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KX 

FIG. 5: The potential difference $0 between the middle of 
FIG. 3: The charge density profile in the slab for several the slab and one electrode as a function of the width W of 
values of the width W at fixed k. the slab. 



^fe" * <x) 




K X 



FIG. 4: The electric potential profile $(2;) for different values 
of the width W of the slab at fixed k. From top to bottom 
kW = 10, 5, 1, 0.16. 



the slab (x = 0) which can be obtained from the previ- 
ous expression by replacing x by W/2 (the term (|4.22c[) 
in the previous equation has the limit — (ln3)/4 when 
x = W/2). Fig. shows a plot of the potential dif- 
ference between the middle of the slab and a plate 
$0 = $(0) - §{±W/2) = $(0) as a function of W . It 
is interesting to know the limit when W — > 00. From 
Eq. we get 



$n = 



W^ao e z K 



7 



(4.23) 



For an asymmetric electrolyte a non-zero potential dif- 
ference between the middle of the electrolyte and any 
plate builds up although both plates are grounded. The 
sign of this potential difference is given by the parame- 
ter X) q Cq9q- This potential difference is a monotonous 
function (increasing or decreasing depending on the sign 
°f laCa) of the width W with an extremum value for 
W — > 00 given by Eq. I|4.23[) . 



It is interesting to comment a few points on the case 
when only one electrode is present, which has been previ- 
ously studied by Aqua using diagrammatic methods. 
In the limit W -> 00, with X = x + W/2, from Eq. (|4"2l|) 
we recover Aqua's expression for the electric potential 



<I>(,Y) -<!>„ 



(4.24) 



Ei(-A) - e A Ei(-3X) 



We can notice that far away from the electrode the po- 
tential behaves as 



HX) - $ 



27T/3 



where we defined 



7 

■kX 



In 3 



1 e 



-si 



*eff = 



27T/3 



E^7 



In 3 



(4.25) 



(4.26) 



This result suggest the following interpretation. If we 
were to understand this result using a mean field lin- 
earized Poisson-Boltzmann equation, we can suppose 
that the electrode has an effective potential <I> c ff given 
by Eq. 1)4.260 . The potential of the electrode, which is 
zero in our case, gets additively renormalized by $ e ff by 
the effect of the fluctuations around the mean field. This 
interpretation follows the same philosophy that the one of 
the theory of the renormalized charge in highly charged 
colloids j3^, |3(j, except that in this case the potential 
renormalization is due to the effect of the correlations 
and not to the non-linear effects of the mean field the- 
ory. If the electrode was at a fixed potential V, the ef- 
fective potential as seen far from the electrode would be 

v + $ eS [3. 

In the spirit of this interpretation, notice that the 
renormalization of the potential $ c ff is positive if 
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J2 a 9aCa is negative, and it is negative otherwise. This 
potential renormalization only occurs for asymmetric 
electrolytes. 

It is interesting to mention that a similar situation oc- 
curs in the charge renormalization of colloids due to the 
non-linear effects in the mean field approach for asym- 
metric electrolytes, although in the other direction. In- 
deed if the charge, say positive, of a colloid is high enough 
to be in a non-linear regime, but small enough to be in a 
non-saturation regime it has been found that the first de- 
viation (quadratic correction) of the effective char ge f rom 
the bare charge have the sign of J2 a q^Ca [37H38ll39j . In 
particular in an intermediate regime the effective charge 
of the colloid could be higher than the bare charge if 
J2 a QaCa has the same sign as the bare charge. This 
analogy is however only qualitative. The context in the 
case of colloids is somehow different from the one consid- 
ered here. In the colloids the renormalization is due to 
the non-linear effects of the mean field approach and here 
we considered the corrections due to the correlations. 



terior of the system. The potential inside the electrolyte 
has the same sign that J2 a q^Ca- 

As possible perspectives to this work let us mention 
the following. Here we studied the first fluctuations cor- 
rections to the mean field description of the electrolyte. 
We considered the most simple situation at the mean 
field level: both plates are grounded and the mean field 
potential is zero everywhere. We choose to study this 
situation in order to show clearly the effects of fluctua- 
tions. However our method can actually be extended to 
study more general situations, for example when a poten- 
tial difference is imposed between the electrodes. In this 
case the mean field description is the Poisson-Boltzmann 
theory studied by Gouy and Chapman Q which can 
be a linear theory if the potential difference is small or a 
non-linear theory if the potential difference is high. On 
the top of this mean field description, a generalization of 
our method can be used to find the first corrections due 
to the microions correlations (whether the mean field is 
linear or non-linear). 



V. SUMMARY AND CONCLUSION 

We have obtained the first corrections due to fluctua- 
tions to the mean field description of an electrolyte con- 
fined in a metallic slab of width W . We considered two 
models to describe the metallic electrodes. The ideal 
conductor model is more tractable, but it neglects the 
fluctuations of the potential inside the electrodes. For 
this (academic) model the disjoining pressure of the sys- 
tem is always positive and it increases if the separation 
W decreases indicating a repulsive force on the metal- 
lic plates by the electrolyte and a stable system. Also 
we confirmed [3(j that for large separations W the dis- 
joining pressure has an algebraic decay in W~ 3 , pd ~ 
[fc s TC(3)/(87r)] W~ 3 . This large-W algebraic finite-size 
correction is universal: it does not depend on the micro- 
scopic constitution of the system. 

For the more realistic model of the good conductor 
electrodes, the behavior of the pressure is completely 
different. There is no algebraic decay in l/W 3 in the 
pressure. Its decay is exponential e~ 2nW at large dis- 
tances. Furthermore the disjoining pressure is negative, 
thus suggesting that there is an attractive force between 
the electrodes. 

We obtained some results for the density profiles and 
the electric potential which are independent of the model 
used to describe the electrodes. We retrieved [llllrillill 
a very interesting behavior if the electrolyte is asymmet- 
ric, in particular if ^ a 7^ 0. In this case the system 
is not locally neutral, there is a local charge density with 
the same sign that J2 a q 3 ( a near the electrodes. Sim- 
ilarly the electric potential is not zero inside the elec- 
trolyte although both plates are grounded: a potential 
difference builds up between each electrode and the in- 
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APPENDIX A: GENERAL EXPRESSION FOR 
THE DENSITIES 

The density can be computed from the grand poten- 
tial using Eq. Q4.1J1 . However to perform the functional 
derivative for arbitrary fugacities Ca( r ) we should find a 
more general expression for the grand potential than the 
one given by Eq. 12.11fl which is restricted to constant fu- 
gacities satisfying the pseudo-neutrality condition (|2.1|l . 
Similar calculations to the one presented here can also 
be found in Refs. pBl liol l4l| in the case of unconfined 
systems. 

In general, the sine- Gordon transformation allows to 
write the gran d partition function without any approxi- 
mation as |l4 124| 

E=^-Jv<f > exp[-S($] (Al) 

with Zg given by Eq. (|2.9|) and the action S given by 
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S(cb) = - 



f^(r)A^r) 



dr. 



(A2) 



Let us define the Gaussian average 

( ...) G = J_ /^(...) e -§/^)[-^M*. (A3) 

Notice that the covariance of the preceding functional 
Gaussian measure is (0(r)0(r')) g = /3 _1 i>(r, r'). There- 
fore the last term of Eq. (|A2|) is very similar to a normal 
ordering, since by definition 



: exp (-i/3q a <j>(r)) 



e 0qiv(r,r)/2 e -if3q a 4>(r) _ ( A4 ) 



However the important difference is that in Eq. (|A2(I we 
subtract the self-energy v°(r, r) for an unconfined system 
not the self-energy v(r,r) for a confined system. As pre- 
viously mentioned this has very important physical con- 
sequences for confined systems. To proceed it is natural 
to define a pseudo-normal ordering as 

:: exp (-i/3q a 4>(r)) :: = e (3?>o(^)/2 e -ifc*(r) ( A5 ) 

and write down the action as 

8e 



S = - 



8tt 



(r)A#r)+£C«(r) 



-i/3q a 4>(r) 



dr. 



(A6) 



As we mentioned earlier in Sec.[n]if we use the Coulomb 
potential v° defined by Eq. (|2.4[1 the self-energy v°(r,r) 



is infinite. In principle one should introduce a short- 
distance cutoff. Without this cutoff both quantities 
e -i/3q a ct>(r) an( j v 0(j^ r ^ m (|A5I> are not properly de- 
fined when taken separately. However their combination 
in the definition (| A5|) of the pseudo-normal ordered ex- 
ponential :: exp (—i0q a (f>(r)) :: is well defined even in the 
case of a vanishing short-distance cutoff. In particular 
later on we will proceed to do an expansion of this expo- 
nential in powers of :: (— i(3q a (j>(r)) n :: . These quantities 
are well defined for a vanishing short-distance cutoff and 
they are of order (T a ) n in the coulombic coupling con- 
stant. 

For arbitrary position dependent fugacities the station- 
ary point of the action S is 4> — —iip with ip solution of 
the mean field Poisson-Boltzmann equation 



4tt 



AV>(r) + — J2( a (r)q a t 



-/3<J„V(r) 



(A7) 



Notice that if one takes constant fugacities satisfying the 
pseudo-neutrality condition ^ Q Ca<Za = the solution 
of the mean field Poisson-Boltzmann equation is simply 
ip(r) = 0. 

Let us return to the general case of position dependent 
fugacities. Expanding the action to the quadratic order 
in (j) around the stationary point leads to S(—iip + </>) = 



S, 



mf - 



Si 



2 ), with 



S m f = S(-iip) 



Be 
8^ 



V-(r)A^(r)-^Ca(r) 



(Zi- 



ffle action evaluated at the mean field solution and 



-0s 
< 

47T 



(r)A0(r)+^(/3 gQ ) 2 C a (r)e 

CX 

I 



-/3q a if>(r 



dr. 



(A8) 



(A9) 



We can now compute the functional derivative 14.1fl with In particular from Eq. (|A7|) we have 
respect to the fugacities to find 



n a (r) = - 



SS, 



mi 



r m , SSl a, 

J <JCa(r) 



(JCa(r) fV<t>e- s i 



(A10) 



S( a (r>) 



4tt 



qj{r-r') (All) 



However we should take special of the terms that depend 
on the mean field ?/>(r) since the latter is a function of 



where 



is evaluated for constant fugacities satis- 



the fugacities via the Poisson-Boltzmann equation i|A7|l . fying the pseudo- neutrality condition (|2.1|) and ip(r) = 0. 
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Then we can write 



9a"DH(r,r') 



(A12) 



with wdh ( r i r') the Debye-Hiickel potential satisfying the 
Debye-Hiickel equation l|4.3|l and the imposed bound- 
ary conditions. Taking this into account we find the re- 



quired functional derivatives evaluated at constant fugac- 
ities satisfying Eq. (|2.1[) and ip(r) = 0, 



and 



SSmf 



SUr) 



(A13) 



SSi 



(r)2;; _^J2q% J '« DH (r / , r):: 0(r) 2 :: dr. 



(A14) 



r 



For constant fugacities the action Si reduces to 
1 f ~(3e 



ai„ = i/ 



4tt 



r)A0(r)+^(/3 97 ) 2 C 7 :: 0(r) 2 :: dr. 

(A15) 



If we define the average 



(...) 



DH 



/£></>(• ..) e - Sl lo 



(A16) 



(A17) 



we have 

/?(::0(r) 2 ::> DH =«DH(r,r)-« o (r,r). 

Then replacing (|A13|) and (|A14[) into Eq. I|A10I) and us- 
ing (|A17I) gives Eq. (|4.2[> for the densities. 



APPENDIX B: ON THE IMPORTANCE OF THE 
BOUNDARY CONDITIONS 



As it was mentioned in section ITTT1 there is an impor- 
tance difference on the behavior of the disjoining pressure 
in the case with ideal conductor boundaries and the case 
when the electrolyte is confined by a walls made of ma- 
terial with a dielectric constant e w [3l]]. Namely, for the 
latter case the pressure has an exponential decay e~ 2KW 
at large separations W whereas in the case presented here 
we found an algebraic decay in 1/W 3 which is further- 
more universal, i.e. independent of the microscopic detail: 
notice that the coefficient of W~ 3 in Eq. I|3.8fl is just a 
number independent of the Debye length k^ 1 and of the 
other microscopic parameters. 

Although the ideal boundary conditions case consid- 
ered here is formally obtained when £^, = 00, this limit 
has a very different behavior than in the case < e w < 00 
(the ideal dielectric boundaries case, e w — 0, is also spe- 
cial, in that case there is also an algebraic decay at large 
separations W~ 3 in the pressure [IJ similar to the one 
found here). This difference is not only present on the 
behavior of the pressure but also in the correlation func- 
tions. For instance, it is known that for dielectric bound- 
aries with < s w < 00 the charge correlation along the 



walls have an algebraic decay as |r 1 |~ 3 where rj_ is the 
direction parallel to the walls |43J, |44j whereas for ideal 
conductor boundaries e w = 00 (or ideal dielectric e w = 0) 
this decay is faster that any power law (see the review 
Ref. and references cited therein). 

In this appendix we show how both kinds of boundary 
conditions can be related and understand the presence 
of the universal term in 1/W 3 in the Dirichlet boundary 
conditions case and its absence in the case of insulating 
boundaries. The following analysis relies on a macro- 
scopic description of the electrolyte in terms of collective 
modes which actually disregards the microscopic detail of 
the system but can give a correct description of some uni- 
versal properties of the system for instance the presence 
of the 1/W 3 term in the pressure. 

Let us consider an electrolyte confined in the slab do- 
main D with separation W and with boundaries made of 
a dielectric material with dielectric constant < e w < 00. 
For simplicity and without lost of generality (divide e w 
by e) we will take the dielectric constant of the solvent 
e = 1, It is well-known that the electric potential $ 
of a linear collective mode of oscillation with frequency 
u> of this electrolyte can be described a Laplace type of 
equation 



X„A$ = 



(Bl) 



inside the domain D, with an effective dielectric constant 
Xu = 1 — uit I LQ 1 and ui p is the plasma frequency (see for 
instance [46|). Outside the domain where the electrolyte 
is confined $ satisfies a Laplace equation A$ = since 
there are no real charges outside. The potential should 
also satisfy the boundary conditions: $;„ = $ out and 
^ui^n^out = X^dn&in, where (9„$ denotes the component 
of (minus) the electric field normal to the boundary. 

One can distinguish between two type of modes. If 
Xui 7^ 0, $ has a vanishing Laplacian inside the domain: 
there are no charges inside the domain. Only at the 
boundaries there are some surface charges. These are 
the surface modes. They represent a system of two par- 
allel walls in the vacuum with fluctuating surface charges. 
Then the contribution of the surface modes to the pres- 
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sure is the same as the one of the Lifshitz theory |3J] 
in the classical (i.e. non-quantum) limit [47L |4S| . This 
contribution comes from the well-known Casimir forces, 
it is attractive and for large separations it is given by 
-k B T((3)/(8TrW 3 ). 

If uj = Up, = and charges inside the domain can 
exist since A$ 7^ inside D is acceptable. These are the 
volume modes which oscillate all at the same frequency 
uj = uj p , the plasma frequency. Due to the boundary con- 
ditions, the potential for volume modes satisfy Neumann 
boundary conditions outside the domain 9„$ out = 0, and 
since $ is harmonic outside and vanishes at infinity, this 
implies that $ = everywhere outside the domain D. In 
conclusion the volume modes represent a system of vol- 
ume charges that satisfy Dirichlet boundary conditions 
for the potential. The volume modes of the electrolyte 
confined by dielectric boundaries are very similar to the 
system we studied here with ideal conductor boundaries. 
If one computes the contribution to the pressure coming 
from these volume modes one will find the same result 



that the one we have found in section ITlI Al namely that 
the pressure has an universal algebraic decay for large 
separations W given by +kBT((3) / (8nW 3 ) . This contri- 
bution gives a repulsive force and is exactly the opposite 
as the one coming from the surface modes. 

Adding both contributions from the surface and the 
volume modes one finds that the algebraic contributions 
to the pressure cancel each other. In conclusion for the 
system confined by dielectric boundaries there is no term 
in 1/W 3 in the pressure as previously noted jjjj. How- 
ever for ideal conductor boundary conditions, a repulsive 
term in 1/W 3 for the pressure is present. 

The above analysis is somehow similar to the one done 
in the comparison between ideal conductor and good con- 
ductor model, in the sense that the absence of the 1/W 3 
algebraic term in the pressure in the dielectric boundaries 
case and the good conductor case is due to a cancellation 
between the Lifshitz term and the one from the ideal con- 
ductor model. However the specific details arc different. 
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